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Lattice distortions in which the axes are locally rotated 
provide an intrinsic source of frustration in anisotropic su- 
perconductors. A general framework to study this effect is 
presented. The influence of lattice defects and phonons in d 
and s + d layered superconductors is studied. 
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I. INTRODUCTION. 

In anisotropic superconductors, local rotations of the 
lattice modulate the order parameter. The most strik- 
ing manifestation of this effect is the influence of grain 
boundaries in high-Tc superconductors |l]jic|,|3|| . Less pro- 
nounced effects are also to be expected from lattice de- 
fects such as dislocations, as they twist the lattice axes 
in their vicinity. Transversely polarized phonons can, in 
principle, also couple to the superconducting order pa- 
rameter through the local rotations that they induce. 

If, along a given direction, there are strains which 
change the orientation of the lattice axes, the equilib- 
rium order parameter must follow that distortion. For 
instance, in a two dimensional square lattice, a param- 
eter with dx2_y2 symmetry, A^2_y2, must change into 
cos(20)A2,2_j,2 — sm{26)Axy after a lattice rotation by 
an angle 6. We express the local order parameter in 
terms of its components in a fixed frame of reference 
which is independent of the orientation of the lattice. 
We expand the local order parameter in terms of the 
functions (A,j.2 _y2 , A^y) , defined in this external frame. 
In principle, the symmetry of the square lattice allows 
for mixing of a pure Ax2_y2 — Aq cos(20) function with 
higher spherical harmonics. Under a rotation, each of 
these harmonics behaves differently, and requires a spe- 
cific description in the external frame of reference men- 
tioned earlier. For simplicity, we will first ignore these 
higher order components. Their inclusion will be dis- 
cussed later, as well as extensions to more complex order 
parameters. 

II. FORMALISM. 

We will now set up an effective Ginzburg Landau de- 
scription of the situation described above. Let us first 
consider a model superconductor with very small coher- 
ence length. At each position, the order parameter will 



match perfectly the orientation of the lattice axes. There 
is no cost in free energy due to misorientations in the or- 
der parameter. We can ignore gradient terms, and use 
an expression which contains only quartic and quadratic 
terms, provided that we insert in them the value of order 
parameter with respect to the local axes. 

The coupling between neighboring regions is described 
by gradient terms. There is a free energy cost associated 
to the inhomegeneities of the order parameter. In the 
present case, the degree of inhomegeneity must be re- 
ferred to the fixed, external frame of reference. There is 
an increase in the free energy of the system when differ- 
ent areas of the lattice are misoriented, even if the local 
order parameter is aligned with the respective axes. A 
particular example of this situation is a grain boundary, 
as described in the introduction. 

The difference between the two ways of describing the 
order parameter, with respect to an external frame or 
to the local one, is determined by the degree of ro- 
tation of the lattice from one region to another. Let 
us define the rotation which takes the axes from one 
point to another by the angle 0. As mentioned previ- 
ously, an order parameter which follows the lattice will 
look, in the external frame, like A^2_y2 in one point and 
cos(2^^)A2,2_J,2 -I- sm{29)Axy in the other. The transfor- 
mation which takes one expression into the other can 
be accomplished by modifying the two component vec- 
tor (A^2_y2, Axy) by means of an operator of the type 

J Adr ^ such that: 

I j Aydr = 26(^\l^ (1) 

and: 

^{A,dr_( cos(26i) sin(26l) \ 
^ \-&m{2e) cos(20) ) ^ ' 

The integral is to be taken along the path in real space 
along which the lattice is rotated. As 9 is the total ro- 
tation along the path considered, we can write, for the 
example considered here: 

Ag = 2(VR)(Ty = AOy (3) 

where V is the gradient operator in real space, R specifies 
the local rotation of the lattice axes, and Oj, is a Pauli 
matrix. The rotation R can be written as S : 

_ dxUy — dyUx 



1 



where u(r) are the local deviations of the lattice node at 
f from equilibrium. Eq. is valid for small deviations 
(strains) . 

The expressions above allow us to relate the changes 
in the order parameter as seen in the two frames, and the 
rotations of the lattice, in terms of the local strains. The 
gradient term in the Ginzburg Landau description, which 
should be expressed with respect to the external frame, 
looks, in terms of the order parameter with respect to 
the local axes: 



The previous analysis can be generalized to three di- 
mensional superconductors with arbitrary order param- 
eters. Local rotations of the lattice can be expressed in 
terms of the parameter: 



Ri 



(6) 



J" 



i=x,y 

aoiT-Tf~y") 



J ^ |(a.-A,,)A, 

i—x,y 



where is the fully antisymmetric tensor. In order 
to determine the transformations induced by these ro- 
tations in the order parameter, we must first decom- 
pose it into spherical harmonics with well defined trans- 
formation properties. Then, a gauge field which per- 
forms rotations in order parameter space is added to 
each component. Finally, as in (|), aU ter ms allowed 
by the lattice symmetry are included in the free en- 
ergy. The procedure is cumbersome, but straightfor- 
ward. For instance, we can add to higher harmonics 
which transform like cos[(4n -I- 2)0],sin[(4n -I- 2)9]. The 
gauge field associated to these components is of the type 
l^xyl"^ + quartic terms (5) Ag = 4(2n -I- l)a^\/{da;Uy - dyU^), where cr^ transforms 

one of these parameters into the other. The gauge field 
induced by the higher harmonics of the order parameter 
are also proportional to the gradient of the lattice rota- 
tions, but with a larger constant of proportionality. The 
term with |Agp in the free energy induces a reduction of 
these components of the order parameter. This reduction 
will be stronger for higher harmonics than for lower ones. 



\^x^-yA 



where J- is the free energy. Tc is a number with di- 
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mensions of temperature, and of order ~^ . Finally, 
flo are proportionality constants which 



play no role in the analysis presented below. We take the 
critical temperatures of the Ax2_y2,Axy order parame- 
ters as different, as appropiate for a square lattice. Equa- 
tion (|^) and the definitions (^) and ^ suffice to study the 
phenomenology of layered superconductors whose order 
parameter is well approximated by [cos(26'), sin(20)]. 

It is easy to see that, when A2.2_y2 and A^y have the 
same critical temperature, and V x A = 0, the effects 
of the field can be made to vanish by performing a lo- 
cal rotation of the order parameter [Q. In this case, the 
system has isotropic superconducting properties. Hence, 
rotations of the underlying crystal lattice leave the (de- 
generate) order parameter unaffected. 

The formalism described here can be expressed in 
terms of a connexion, which depends on the lattice orien- 
tation. The definition of parallel transport of a vector, in 
our case (Aa,2_y2, Aj;j,), needs to be modified by the ro- 
tation of the lattice. The usual derivative is changed into 
a covariant derivative |^ . This technique lies at the basis 
of the extensive work done in topological phases ||^ , and 
has found use in many fields in condensed matter physics 

In addition to the effects described by {m, there may 
be an explicit coupling of the order parameter to the 
lattice strains. The simplest coupling is proportional to 
J J2i Ciil Aa;2_j,2 p, where J2i defines the local com- 
pression, or expansion, in terms of the symmetric part of 
the strain tensor. This term can also exist in a s wave 
superonductor. A coupling of this type is not negligible 
in the high-Tc compounds, due to the strong dependence 
of Tc on pressure 0] . 



III. STATIC DISTORTIONS. 

We now consider specific lattice distortions within the 
layers. The simplest static defect is a dislocation. The 
strains in the vicinity of a dislocation can be calculated 
from the theory of elasticity. If the Burgers vector is b, 
then g. 



dx 



dyUx 



br 



(7) 



where we are neglecting a possible anisotropy in the elas- 
tic constants of the lattice. We now assume that the A^y 
component is negligible, which is justified if, for instance, 
the parameter T^*' in (^ is less than zero. Then, the cou- 
pling of the lattice rotations to A^2_y2 is due to terms 
quadratic in the gauge potential, and we obtain: 



O'xyTc^Q 



|A^2_,, 



(8) 



and lAI 



The effect of this term is to reduce 



|6P 

the effective critical temperature near the dislocation. It 
suppresses a possible local enhancement of Tc due to the 
explicit dependence of the critical temperature on pres- 
sure 1^ . In the presence of a magnetic field, the normal 
cores of vortices will be attracted towards the disloca- 
tions, to take adventage of the supression of Tc, leading 
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to pinning. Near the core of the dislocation, r '-^ ^q, we 
can estimate the reduction of Tc (see below) as ^ . 
This effect can be significant in a material with a short 
coherence length, such as the copper oxides. 

A grain boundary can be modelled as an array of dis- 
locations. The gauge fields induced by each of them can 
be added, to obtain the effect of the whole boundary. Af- 
ter some algebra, we find that, if the rotation of the axes 
across the boundary is not large: 

^ ^ sin\0) [1 - cos ) cosh {^)]' + sin^ {^) 



2o? 
TT^ sin^(6l) 



[cosh (^)- cos (f)] 



(9) 



where a is the lattice constant, is the angle between the 



axes at both sides of the boundary, and a 
(a and " 



From 



sin(e) ■ 

we can estimate the reduction of the critical 



temperature at the boundary. We assume that a^2_y2 



ao, and fc - T^'"^ . Then, ^ 



where Tc = ■ This calculation is rigorously valid 

when sin(0) <^ 1, that is, when strains change slowly. A 
significant reduction of Tc can be expected near a twin 
boundary, which will act in a way similar to a weak link. 
It would be interesting to check this experimentally [f2| . 



IV. DYNAMIC ROTATIONS. 

The orientation of the lattice can also be modified dy- 
namically, by transversely polarized phonons. A longitu- 
dinal phonon induces local compressions and dilatations 
of the lattice. A transverse phonon does not change the 
volume of the unit cell, but induces a shear deformation 
which changes from position to position. This deforma- 
tion can be viewed as a local rotation of the axes. In 
a material with a complex unit cell, there are also opti- 
cal phonons which rotate locally the axes. The rotations 
can be analyzed by means of the formalism discussed be- 
fore. The only difference with a static distortion is that 
the gauge field which accounts for the lattice rotations 
depends on dynamic displacements. The time scales as- 
sociated to the superconducting response are of the or- 
der of the inverse plasma frequency. This frequency is 
much higher than typical phonon frequencies, so that the 
changes in Tc induced by phonons can be studied as if 
the lattice distortions were static. 

Transversely polarized acoustic phonons give rise to 
strains proportional to the gradient of their amplitudes. 
Hence, we expect that the Fourier components of R are 
given by: 



k^ X Ufc 



(10) 



where is the amplitude of a phonon of momentum k, 
and k±_ is the component of k perpendicular to Ufc. For 



transversely polarized phonons, k is always perpendicular 
to Ufc, and we do not need to take this restriction into 
account. 

The gauge potential in (||) has a finite average, related 
to the thermal (or quantum) average value of the phonon 
amplitudes. It is given by: 



|A,p^^(P|(kxufc)p} 



(11) 



We are analyzing materials which show strong dif- 
ferences in the superconducting properties along the in 
plane and out of plane directions. Phonons, however, 
need not be so anisotropic. We assume that they are 
described by an isotropic sound velocity, c, and a Debye 
frequency ojd ~ — , where a is the lattice constant, which 
we take to be of the same order of magnitude in all di- 
rections. The three dimensional integral implicit in ( pT| ) 
is extended to those modes such that Tick < ksT and 
a < k^^ and < k^^. Note that the Ginzburg-Landau 
des crip tion used here is valid at length scales greater than 
^0 ||ll|l • For simplicity, we take the Bose occupancy of a 
mode of frequency uj to be None of these approx- 

imations will alter the order of magnitude of the effect. 
We can distinguish three regimes: i) if fcsT <C hcS,Q^, the 
upper cutoff in (|ll| ) is due to the Bose factor, ii) when 
7ic^(7^ <C ksT <^ fica~^, the upper cutoff is fc|| ~ 
and kz ^ and iii) if fica~^ <C ksT, the maximum 
value of kz is a~^. We can now estimate the reduction 
in the (average) value of the critical temperature in the 
same way as in the previous section. We find: 



ATc 
Tc 



h kBTc 



fceTc « 



hup a 
Co 



(12) 



where p is the density. This reduction of the critical tem- 
perature should not be negligible in the high-Tc com- 
pounds, as ^0 is not very different from a, the lattice 
constant. 

In addition to a reduction in the average value of Tc, 
we can also analyze the fluctuations in this value. We are 
treating phonons as a source of static disorder, which re- 
duces Tc. The value of A^ fluctuates, changing the local 
values of the critical temperature. If these fluctuations 
take place on scales larger than ^q, a distribution of Tc's 
will be observed. As the changes in Tc go as: 



ATc^A2^^|A:|V|fifeP 



we flnd that: 



(ATc(f)ATc(f'))^^k« 



zk(r-r')|^j^|4 



(13) 



(14) 



and: 
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(ATc 



T, 



7 hujjja 

?0 



(15) 



by comparing ( p^ and this expression, one can see that 
fluctuations cannot be neglected. Finally, the typical 
length scale for fluctuations inside the planes is 



tie 



w 



^0 (16) 



where W is the electronic bandwidth. Hence, we ex- 
pect large fluctuations in small scale domains. Transport 
measurements will be unable to resolve this structure, 
and the average Tc will be observed. 

Transverse optical phonons of frequency uJopt- give an 
additional reduction of Tr, 



AT, 



a- 



(17) 



opt 



where {v?) denotes the mean square displacement of the 
phonon. It is, typically, a fraction, 10"^ — 10^^ of a'^. As 
in the previous case, we assume that the superconducting 
order parameter has a short coherence length in the out 
of plane direction. There is a contribution like (jlj) for 
each optical phonon transversely polarized in the plane. 

The frustration of the order parameter in anisotropic 
superconductors by dynamical vibrations implies that T^. 
cannot be much larger than the frequency of the phonons 
responsible of the effect. 



V. COUPLING BETWEEN DISTORTIONS AND 
CURRENTS. 

The Ginzburg Landau expression (^) contains terms 
which are linear in Ag, which couple one component of 
the order parameter and the spatial derivatives of the 
other (note that Ag is an operator which rotates one 
component into the other). These derivatives are differ- 
ent from zero in the presence of a supercurrent, because 
of the variation in the phase. In the case described by 
(^, however, these terms do not contribute, as lS.xy ~ 0. 

Terms of this type play a role in materials where the 
superconducting properties differ significantly along the 
lattice axes. Let us consider an orthorhombic supercon- 
ductor, where the order parameter is s + rf. In addition 
to the terms in (||) we also have: 



'tho 



2 
2 



dyA^d^Aly + cc. 



I A, 



a,d(r - T'*) 
c.c. H- quartic terms 



(18) 



where As is the s component of the order parameter. 
For the sake of simplicity, we ignore terms induced by 
the orthorhombic symmetry which are quadratic in the 
gradients |l7j. These terms, needed for a complete de- 
scription of the anisotropy in the penetration depth, do 
not modify qualitatively the picture discussed in the pre- 
vious sections. 

In the absence of internal lattice rotations, A^y — 0. 



The first term in (^^ leads to a difference in the penetra- 
tion depth along the two lattice axes. Local rotations of 
the lattice can be incorporated by adding the gauge field 
(^) to the derivatives of the d components of the order 
parameter in (|l^). Then, assuming that the equilibrium 
order parameter has both Aa;2„g2 and A^ components, 
we obtain terms of the type: 



a^AsAgA*2 



dyA.A^A 



^ ~y 



c.c. 



(19) 



This expression gives a contribution if VA^ ^ 0, and, in 
particular, in the presence of a current. 

The changes in a current carrying state induced by 
lattice rotations are best studied within London's frame- 
work, which is valid if the magnitude of the order param- 
eter does not change. We define the London tensor, with 
the same symmetry properties of the perfect lattice, such 
that j = m~-'-Aeiec: where A-^iec is the electromagnetic 
vector potential. For the study of currents within the 
layers, we have: 



(20) 



where Xx , Ag are the penetration depths along the two 
axes, and c is the velocity of light. London's expression 
for the free energy is |13]: 



1 

8^ 



IV X A, 



elec I 



! ^ 

'y 2 



(21) 



The deviation of tetragonal symmetry is given by A^ ^ — 

xf^^^\AsAl._^.\. 

Taking into account that the gauge field related to the 
lattice distortions. A, is itself a derivative, we can inte- 
grate by parts in eq.(|l9|), to obtain: 

ATl = ^ / a.A,R9gA:._^. 



oc (A-^ - Ag ^)A^A^ j jxidxuy - dyUx)jy (22) 

This equation can be interpreted as a rotation of the 
London tensor, (pO|), induced by the distortion of the lat- 
tice. The currents adjust to the directions of the local 
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principal axes. At a grain boundary, for instance, AeZec 
must be continuous. The component of the current par- 
allel to the boundary must be discontinous, leading to a 
kink in the direction of the total current, as the London 
tensor which gives the current at both sides is discon- 
tinuous |l4-16|. If, on the other hand, the boundary 
conditions prevent the current from acquiring a trans- 
verse component, the vector potential shows a kink, and 
a perpendicular magnetic field is generated at the bound- 
ary, as sketched in fig(0 [|l8|. This picture, derived from 
London's theory, is valid at distances comparable to the 
penetration depth. The discontinuity at the boundary is 
rounded off at shorter scales. 




FIG. 1. Currents and electromagnetic vector potentials 
near a twin boundary, in a slab much narrower than the pen- 
etration depth. When the current is constant, and flows par- 
allel to the boundaries, the vector potential must have a dis- 
continuity at the twin. This discontinuity induces a magnetic 
field in the direction perpendicular to the slab. 

The present analysis implies that currents can also be 
deflected by other types of static defects, and by phonons. 
The currents flowing near a dislocation, and the instan- 
taneous current in the presence of phonons, are sketched 
in fig.(||). Magnetic fields will be induced. It is interest- 
ing to note that phonons can modulate the currents and 
give rise to electromagnetic fields. 





a) 



b) 



FIG. 2. Sketch of the a currents flowing near a dislocation 
(a), and instantaneous current in the presence of a transverse 
acoustic phonon (b). 

So far, we have only considered the role of in plane dis- 
tortions in layered superconductors. If we include rota- 
tions around in plane axes, we must allow for components 



like A322_r2 , Aa;^, in the order parameter. Because 
of the strong anisotropy of these compounds, components 
which transform like higher powers of are also to be 
expected. In analogy with (|2^) , these components lead to 
contributions of the type / ji{diUz — dzUi)jz, where the 
direction i lies in the planes. Note that this effect should 
also be present in layered superconductors with s sym- 
metry within the layers. The magnitude of the J 

coupling is proportional to A ,Ja ■ 

I ' II 




a) 




FIG. 3. (a) Sketch of the currents around a vortex pinned 
by a screw dislocation. Dashed line is the dislocation axis, 
(b) Currents flowing along a twisted monocrystaline whisker. 

A particular defect which generates rotations in the 
out of the plane direction is a screw dislocation, which 
seems to be common in sputtered films |19 . In cylindri- 
cal coordinates, we have that only deUz — ^ '^s different 
from zero |8|. Hence, the free energy of the superconduc- 
tor has a term like J jgjz ■ A current flowing around the 
dislocation induces a current parallel to the axis of the 
dislocation. This situation may arise if a vortex is pinned 
by a screw dislocation. A similar effect can be realized in 
a monocrystalline whisker subject to shear strains. Then, 
dzUe 7^ 0, leading again to a J jejz coupling. A current 
flowing along the axis of the whisker induces a tangential 
component, giving rise to a magnetic field in its inte- 
rior. A schematic picture of both situations is sketched 
in fig.®. 



VI. CONCLUSIONS. 

We have presented a formalism to study the infiuence 
of lattice rotations on anisotropic superconductors. The 
analysis has been particularized to d-wave layered sys- 
tems, but it is general enough to be applicable to more 
complex lattices and order parameters. We expect the 
scheme outlined here to be also useful for the study 
of heavy fermions and other materials which show an 
anisotropic order parameter. 

For the systems considered here, we find that 
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i) Lattice distortions couple directly to an anisotropic 
order parameter, irrespective of other mechanisms, like 
the dependence of Tc on pressure. 

ii) Higher harmonics of the order parameter are more 
strongly frustrated by local rotations Q . 

iii) Lattice defects, such as dislocations, reduce locally 
Tc- In the presence of a magnetic field, the vortex cores 
will be attracted to these regions, giving rise to pinning. 

iv) Transversely polarized phonons also reduce Tc. 
The value of the critical temperature cannot be much 
larger than the typical frequency of these phonons. 

v) Local rotations deflect supercurrents in supercon- 
ductors with inequivalent lattice axes. In particular, the 
current through a grain boundary is deflected. 
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